Kinematic dynamo simulations of von Karman flows: application to the VKS 

experiment 



A. PinterQ B. Dubrulle, and F. Daviaud 

Service de Physique de VEtat Condense, CEA, IRAMIS, 
CNRS/URA 2464, CEA-Saclay, F-91911 Gif-sur-Yvette, France 
(Dated: Eur. Phys. J. B 74, 165-176 (2010)) 

The VKS experiment has evidenced dynamo action in a highly turbulent liquid sodium von 
Karman flow [R. Monchaux et al., Phys. Rev. Lett. 98, 044502 (2007)]. However, the existence and 
the onset of a dynamo happen to depend on the exact experimental configuration. By performing 
kinematic dynamo simulations on real flows, we study their influence on dynamo action, in particular 
the sense of rotation and the presence of an annulus in the shear layer plane. The 3 components of 
the mean velocity fields are measured in a water prototype for different VKS configurations through 
Stereoscopic Particle Imaging Velocimetry. Experimental data are then processed in order to use 
them in a periodic cylindrical kinematic code. Even if the kinematic predicted mode appears to be 
different from the experimental saturated one, the results concerning the existence of a dynamo and 
the thresholds are in qualitative agreement, showing the importance of the flow characteristics. 

PACS numbers: 47.65.-d, 91.25. Cw 



I. INTRODUCTION 

Dynamo action is an instability by which mechanical 
energy is transformed into magnetic energy for R m > R^ 
where R m is the magnetic Reynolds number. Since Lar- 
mor's prediction concerning the solar dynamopQ, a lot 
of work has been devoted to this field, including both 
theories and numerical simulations. In the last decade, 
the interest for dynamo has been renewed by the posi- 
tive results obtained by the Riga [2] and Karlsruhe [3] 
fluid dynamo experiments. These experiments rely on 
analytical predictions by Ponomarenko [1] and Robert 
[5] and have been designed accordingly. In particular, 
the generated turbulent sodium flows are organized so 
that the instantaneous flow remains always close to the 
mean flow, allowing the kinematic dynamo predictions 
of the thresholds and of the neutral modes (respectively 
oscillatory and stationary) to be very accurate. 

A new step was passed very recently in the von Karman 
Sodium (VKS) experiment. Different dynamos have 
been observed inside a flow of liquid sodium at large 
Reynolds number (Re ~ 10 7 ) [5HS]. Depending on the 
the control paramaters, statistically stationary dynamos 
or dynamical dynamo regimes including field reversals 
have been reported. Under these conditions, the flow is 
fully turbulent, with many large and small scale fluctu- 
ations, and the instantaneous flow is very different from 
the mean flow [101 E]. The dynamo instability taking 
place in the VKS experiment is thus at variance with 
the classical instabilities taking place in laminar flows or 
in turbulent flows remaining very close to their mean. 
Nevertheless, many kinematic dynamo simulations have 
been performed to optimize this flow, using the synthetic 
MND flow [TJ] or time-averaged flows measured in water- 
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prototype experiments [TH |5D] . 

Different observations have been made that can help 
elucidate the exact mechanism by which the dynamo in- 
stability develops in the VKS experiment: i) in the ex- 
actly counter-rotating case, dynamo action is observed 
above a critical magnetic Reynolds number Rf n 32 and 
the dominant saturated mode is an axisymmetric m=0 
mode in cylindrical coordinates [6]; ii) so far, dynamo has 
only been observed with iron propellers [3] ; hi) in a given 
experimental set up with iron disks, dynamo has been 
observed when counter-rotating disks in (+) 'non scoop- 
ing' direction (with respect to the propellers blades), but 
no dynamo has been observed when counter-rotating the 
disks in the (-) 'scooping' direction [3]. 

Observation i) shows that the dynamo mechanism is 
not the simple linear response to the mean flow alone, 
since the latter is axisymmetric and produces a trans- 
verse m = 1 mode (m denotes the azimuthal wavenum- 
ber) above R c m ~ 43 as indicated by kinematic dynamo 
simulations [14] • It has been suggested in [15] that the 
m=0 mode of the VKS experiment is in fact generated by 
non-axisymmetric fluctuations through an alpha-omega 
mechanism taking place within the blades of the pro- 
pellers. This suggestion has been confirmed by recent nu- 
merical simulations |16| showing that the addition of an 
(empirical and large) alpha term in the induction equa- 
tion changes the dominant m=l mode into a dominant 
m=0 mode with a dynamo threshold close to the observed 
experimental critical value R c m ~ 32. Another interest- 
ing suggestion has been made by |17| . namely that the 
m=0 mode could arise from a secondary bifurcation of 
the m=l mode, since the latter breaks the original ax- 
isymmetry of the velocity field. 

Observation ii) shows that magnetic boundary condi- 
tions have a large impact on the dynamo instability. High 
magnetic permeability conditions are known to decrease 
the dynamo threshold for the Roberts and Ponomarenko 
flows |18j and more recently for von Karman type flows 
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[!§] . For iron disk, it is suggested that the dynamo 
threshold could be reduced by a factor as large as 2 with 
respect to stainless disks. Since it has not been yet possi- 
ble to get beyond R m = 50 in the VKS experiment, this 
could explain why dynamo has not been observed with 
stainless steel disks, while it is observed around R m = 32 
with iron disks. Other kinematic simulations |20| have 
shown that the sodium flow behind the disks can lead to 
an increase from 12% to 150%. Using iron disks can thus 
also screen magnetic effects of the flow behind the disks. 
PARLER des DISQUES CAMEMBERTS? 

Observation iii) shows that the dynamo is not pro- 
duced only by the rotation of magnetized iron disks in 
a conducting fluid but depends on the flow characteris- 
tics: it is a real fluid dynamo. However, not much has 
been done so far to try and exploit this observation. In 
particular, no kinematic dynamo simulations have been 
performed so far with experimental flows including an 
annulus as in the VKS experiment or with propellers 
counter-rotating in the (-) 'scooping' direction. 

In the present paper, we use Stereoscopic Particle 
Imaging Velocimetry (SPIV) measurements obtained in 
von Karman water (VKE) experiments modeling the 
VKS experiment to perform kinematic dynamo simula- 
tions. We investigate the role of the mean field structure 
on the dynamo threshold, with special attention to the 
rotation sense. Because of the axisymmetry, the exhib- 
ited kinematic mode is, as expected, different from the 
experimental saturated one. However, we show that the 
kinematic dynamo thresholds qualitatively vary like the 
observed dynamo thresholds and discuss how these re- 
sults can be useful to understand the actual mechanisms 
at the origin of VKS dynamo. 

Our paper is organized as follows: In Sec. II, after 
the problem formulation, we describe the experimental 
measurements and the numerical method. In Sec. Ill, 
we focus on the processing procedure of the measured 
velocities and perform a comparison with previous esti- 
mates based on Laser Doppler Velocimetry (LDV) mea- 
surements. Results for dynamo action are presented in 
Sec. IV and compared with the VKS experiments. Fi- 
nally, we discuss the consequences of these results for the 
VKS dynamo. 



II. METHODS 

A. Problem Formulation 

In the present paper, we focus on the kinematic dy- 
namo problem: we consider a given velocity field u, e.g. 
as measured in a water experiment, and check its ability 
to sustain dynamo action by solving the magnetic induc- 
tion equation 

^ = Vx (uxb) + — Ab, (1) 
ot fia 



where b denotes the magnetic field with V • b = 0, fx and 
a are respectively the magnetic permeability and the elec- 
trical conductivity of the considered medium, and A is 
the Laplacian operator. Since u is given, the equation is 
linear in b. There are two possible long-time behaviour 
for the magnetic field: either it is decaying and there is no 
dynamo, or it is increasing, and there is dynamo. Growth 
or decay is in general exponential (although there can be 
transient algebraic growth, due to non-normality of the 
operator, see [21]). To quantify the growth, we introduce 
the growth rate a = d t In (-B 2 ) that is based on the mag- 
netic energy (-B 2 ), where the average is taken over the 
spatial domain. 

Since the equation is linear, no saturation is observed 
in the kinematic approach. This technique is therefore 
unable to detect subcritical or finite amplitude dynamos 
[22] or to study the saturation regime [25H27] as in nu- 
merical flows such as the Taylor-Green flow.CITER SIM- 
ULATIONS FOREST??. However, this approach gives 
an estimation for the dynamo threshold as a function of 
the mean velocity field. 



B. Experimental measurements 

Geometry and experimental setup. — The setup corre- 
sponding to the VKS experiment is described in details 
in [5] and sketched in Figure [lja) . A cylindrical vessel 
with a fixed aspect ratio of L = H/R = 1.8, where R 
(resp. H ) denotes the radius (reps, the length) of the 
vessel which is filled with liquid sodium. We arrange 
the coordinate system such, that the axial extension of 
the cylinder is situated in the interval z E [—0.9,0.9]. 
The fluid is stirred by two independent rotating so called 
TM73 impellors [14] , at frequencies fx and ji. The im- 
pellers are fitted with curved blades (see Fig. [I] (b)). 
Therefore, the flow properties depend on the sense of ro- 
tation, specifically on whether the blades are presenting 
their concave ((-) rotation) or convex ((+) rotation) side 
to the flow during impeller rotation. In the special case 
where /i = /2 > (resp. fx = fi< 0) with the sign con- 
vention of Fig. [lj the impellers fullfill exact (+) (resp. 
(-)) counter- rotation. An optional concentric annulus can 
also be inserted at mid height and has been used in the 
first dynamo runs. This annulus reduces large scale fluc- 
tuations |11) and can also have an impact on the mean 
flow, in particular on the toroidal to poloidal ratio T (see 
below). In the VKS experiment, an additional layer of 
resting sodium is present in the outer part of the experi- 
ment. 

The VKE experiment is the exact half-scale water 
model of the VKS experiment, except that the experi- 
ment is realized vertically (vertical z axis) whereas the 
VKS setup is horizontal. The cylinder radius and height 
of the VKE vessel are R — 100 and 240 mm respectively. 
We have used TM73 impellers of diameter 150 mm fitted 
with radial blades of height 20 mm and curvature radius 
92.5 mm. The inner faces of the disks are H = 180 mm 
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apart. We have worked with different vessel geometries, 
allowing the insertion of an annulus - thickness 5 mm, in- 
ner diameter 170 mm - in the equatorial plane. Note also 
that in the VKE experiment, the layer of resting fluid is 
absent. The velocity field is therefore only measured in 
the flow region, see Fig. [I] (a). 

Velocity measurements are performed using a DAN- 
TEC SPIV system, providing the three components of 
the velocity field on a plane through a 95 x 66 points 
grid [25]. In all that follows, the measurement plane is a 
meridian plane containing the rotation axis. The velocity 
fields are averaged over time series of 5000 regularly sam- 
pled values (sampling frequency between 1 and 15 Hz). In 
the exact counter-rotating regime, the system is axisym- 
metric and symmetric with respect to any rotation R v 
of 7r around any radial axis in its equatorial plane. The 
mean flow is composed of two toroidal cells separated by 
an azimuthal shear layer and two poloidal recirculation 
cells. Figure [2] shows the velocity fields obtained in the 4 
studied configurations: (+) and (-) rotation sense in the 
presence or the absence of an annulus in the equatorial 
plane. Note that in the exact counter-rotating regime, no 
turbulent bifurcation towards a one cell state is observed 
with the TM73 impellers contrary to what is reported 
for impellers with higher curvarture blades in [55 J. The 
bifurcation is observed in the (+) sense for a rotation 
number 9 = (f 1 — /2)/(/i + f2) = 0.1 without annulus 
and 8 = 0.17 with annulus (TT] . 

Control parameters. — Choosing the magnetic field dif- 
fusion time /iai? 2 as the time scale, and the maximum 
absolute value of the flow velocity |u| max as velocity unit, 
the induction equation can be expressed in dimensionless 
form: 



— = R„y x (u x b) + Ab, (2) 
ot 

where R m — /icri?|u| max is the magnetic Reynolds num- 
ber. The critical magnetic Reynolds number R c m for dy- 
namo onset is then defined such that a{R c i nl ) — 0. An- 
other important parameter to characterize the velocity 
field is the poloidal to toroidal ratio: 




(T) = dz \u e (r,z)\rdr, (3c) 

J-L/2 JO 

which is known to have a great impact on the dy- 
namo threshold [13l [14] . Here we denoted the compo- 
nents of the velocity field in cylindrical coordinates as 
u = (u r ,u e ,u z ). 



C. Numerical code 

Because the system is axisymmetric, the time-averaged 
velocity field is supposed to be axisymmetric. We inte- 
grate the induction equation using an axially periodic 
kinematic dynamo code. We only provide here its brief 
description, details can be found in [30]. The code 
is pseudospectral in the axial and azimuthal directions 
while the radial dependence is treated by a high-order 
finite difference scheme. The numerical resolution corre- 
sponds to a grid of 81 points for one wavelength in axial 
direction, 4 points in azimuthal direction (corresponding 
to azimuthal wave numbers m = 0,±1) and 31 points in 
radial direction for the flow domain. The time scheme is 
second-order Adams-Bashforth with a typical time step 
of 10~ 5 . 

Electrical conductivity and magnetic permeability are 
homogeneous and the external medium is insulating. Im- 
plementation of the magnetic boundary conditions for a 
finite cylinder is difficult, due to the nonlocal character of 
the continuity conditions at the boundary of the conduct- 
ing fluid [19) . In contrast, axially periodic boundary con- 
ditions are easily formulated, since the harmonic external 
field then has an analytical expression. Indeed, the nu- 
merical elementary box consists of two mirror-symmetric 
experimental velocity fields in order to avoid strong ve- 
locity discontinuities along the z axis [H] . Note that we 
do not take into account the effet of the flow behind the 
propellers. This effect has been studied in details in other 
kinematic simulations [2D] . 

Finally, we can act on the electromagnetic boundary 
conditions by adding a layer of stationary conductor of 
dimensionless thickness w, corresponding to sodium at 
rest surrounding the flow exactly as in the VKS experi- 
ment in Cadarache [DJ. This extension is made keeping 
the radial resolution of the experimental measured ve- 
locity field constant (31 points in the flow region). The 
velocity field we use as input for the numerical simula- 
tions is thus simply fit in an homogeneous conducting 
cylinder of radius 1 + w: 

u = Ufl ow for < r < 1, (4a) 
u = for 1 < r < 1 + w. (4b) 

Former investigations [T3] showed a decrease of the criti- 
cal magnetic Reynolds number for increasing conducting 
layer size, with a saturation above w = 0.4. Therefore, 
we fixed w = 0.4 (meaning altogehter 43 points in radial 
direction) in order to save computing time. Some cal- 
culations with w — 0.6 (49 points) have also been done, 
when comparing with older computations based on LDV 
measurements. 

Figure [3] gives an example of the total velocity field 
used in kinematic simulations, consisting of an experi- 
mentally measured velocity field Ufl ow and a conducting 
layer with w = 0.6. 
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III. DATA PREPARATION AND TESTS 
A. Preparation of measured velocities 

Due to experimental measurement constraints, it is 
impossible to determine the field in the whole cylinder 
cut (r, z) — ([0, 1], [—0.9, 0.9]): i) in the propeller regions, 
the presence of blades prevents the SPIV measurements; 
ii) close to the cylinder, reflection effects spoil any ve- 
locity measurements; iii) when an annulus is inserted 
at mid-height (z — 0), measurements cannot be per- 
formed in the thin area behind and around the annu- 
lus. Therefore a typical data set covers only the area 
(r,z) as ([0, 0.9], [—0.6, 0.6]). This means, that approxi- 
mately more than a third of the velocity field is unknown 
and has to be determined by boundary conditions and 
spline extrapolation. The radial adaptation has always 
been performed before the axial one. The preparation 
methods of the used SPIV measurements are explained 
in detail in the following. 

Axial Symmetrization. — For the perfect counter- 
rotating setup an axial symmetrization has been per- 
formed with height mid as reflection plane. After this 
procedure u r is strictly symmetric, whereas uq and u z 
are antisymmetric. 

Axial boundary conditions. — The axial velocity field 
should vanish at the blade border, therefore we im- 
posed the condition u z (z = ±0.9) = 0. The two other 
components fullfill a vanishing derivative in axial direc- 
tion at the blades, meaning due/dz(z = ±0.9) = and 
du r /dz(z = ±0.9) = 0. 

Extrapolation in axial direction. — Close to the blades 
the fields are unknown and even not all the data which 
are measured are suitable, meaning that inappropriate 
values have to be removed. This lack is repaired by a 
spline extrapolation. The aforementioned boundary con- 
ditions help us to implement a good fit. The procedure at 
the top and the bottom blade is the same. The number of 
points which have to be removed at a blade are denoted 
by nrp (= "number of removed points") in the following 
text. We compared results for different values of nrp in 
order to find a trustable result with a respective error es- 
timate. Altogether between 11 and 15 points have to be 
extrapolated at a blade corresponding to a nrp between 
2 and 6, meaning that for the first 9 points no measured 
values are available. The number of points used to deter- 
mine the spline function has no significant influence on 
the result, as far as more than 3 are taken into account. 
For security reasons we always took 6 points. 

Interpolation in the annulus region. — When the an- 
nulus is present, there are typically 8 (10 % of the data 
set) measured points around the mid-height that are ab- 
sent or spoiled. They have been removed and replaced 
by spline interpolated values, which is a good natured 
operation. 

Radial boundary conditions. — The radial velocity field 
vanishes at the cylinder border, therefore we imposed 
the condition u r (r = 1) = 0. The two other compo- 



nents satisfy duQ/dr(r = 1) = and du z /dr(r = 1) = 0. 
In the cylinder center the radial and azimuthal compo- 
nents are zero, meaning u r (r — 0) = u@{r = 0) = 0, 
whereas the axial field has a disappearing derivative 
du z /dr(r = 0) = 0. 

Extrapolation in radial direction. — Close to the cylin- 
der the last two points have been replaced and an addi- 
tional point has been added to reach the cylinder border. 
They are the result of a spline extrapolation fullfilling the 
aforementioned boundary conditions using two points for 
the function construction. 

Check with LDV measurements. — LDV measurements 
are free of certain deficiency of the SPIV measurements, 
especially near the cylinder boundary, where reliable ve- 
locity fields can be obtained. It is therefore interesting to 
test the validity of our SPIV field preparation by compar- 
ison with LDV fields obtained previously by Ravelet et al. 
(13) . Figure Q compares the velocity fields of LDV with 
the new SPIV measurements for fixed radial coordinate 
cuts explained in Fig. [3j One cut is close to the vortex 
center, the other one far away. A setup with (+) counter 
rotation has been used without annulus. In addition, the 
processed data, used in the code, with the aid of sym- 
metrization and spline extrapolation for different nrp are 
shown. The difference between the two curves is mainly 
in the azimuthal field «e nearby the blades. Note, that 
the large deviations of the two processed curves for the 
u z field in the regime r ~ [—0.6,0.6] is only due to the 
performed scaling described in the caption. The SPIV 
dataset is finer and overall the comparison is quite satis- 
factory. 

Corresponding cuts for fixed axial coordinate are 
shown in Fig. [5] for planes in the vortex region and in the 
pure inflow layer separating the vortex pair. The agree- 
ment between LDV and SPIV measurements is good even 
if a small deviation can be observed on u r . 

AJOUTER PHRASE CONCLUSION: AVANTAGE 
SPIV/LDV 



B. Influence of the annulus 

In order to study the influence of the annulus, we have 
measured the corresponding time-averaged velocity field 
for different nrp (cf. Fig. [6j. The velocity fields cor- 
responding to the 4 configurations that are used in the 
kinematic dynamo simulations are then shown in Figures 
□ andEl 

In the (+) rotation case, the axial gradient of u r , ug 
and u z are larger in the region of the shear layer. The 
variation with r of the velocity fields looks very similar 
except for u r at mid-radius. However, the poloidal to 
toroidal ratio T, does not reflect this difference: T = 0.8 
in both cases. 

In the (-) rotation case, we were not able to perform 
such a detailed comparison, due to an unexpected ex- 
perimental problem: the locking of large scale vortices 
of the middle shear layer at the position of the small 



5 



cut performed through the annulus to enable SPIV mea- 
surements (cf. t llj for details). This locking does not 
occur in the (+) case, resulting in an axisymmctric time- 
averaged velocity field. In the (-) rotating case, how- 
ever, the locking artificially breaks the axisymmetry of 
the time-averaged field, preventing a priori the compar- 
ison with the case without annulus. Computation of T 
reveals a difference between the left (r = 0.89) part of 
the SPIV plane where the vortices are trapped and the 
right (r = 0.62) part where there is no trapping. Other 
measurements show that only this right part is relevant 
for the VKS experiment in which the annulus is cut-free. 
In the following, we use this right part and symmetrize 
the velocity field accordingly. The corresponding T in 
(-) rotation sense are thus 0.62 with annulus, and 0.49 
without annulus. 

To sec the influence of T on dynamo action for these 
velocity fields, we have also built a synthetic velocity field 
by taking the velocity field measured in the (-) rotation 
case without annulus, and rescaling its poloidal compo- 
nents to reach T = 0.8. This synthetic field is represented 
in Fig. [7]and|8] The variations of its velocity components 
with z and r are very similar to those obtained for the 
(+), contary to the variation of the (-) case with annulus. 



compared also our computations with former LDV results 
presented in [14] . The used fields have already been pre- 
sented in figures |3] and [5] for a (+) rotating setup without 
annulus. Therefore we computed for different nrp values, 
namely 2, 4 and 6, the magnetic energy growth rate as a 
function of R m . The results are presented in Fig. [9] and 
compared with the corresponding LVD curve. The com- 
puted critical magnetic Reynolds number of SPIV varies 
in a range R c m = 41.6 — 47.6, the highest value corre- 
sponding to nrp=6. The comparison of velocity profiles 
in Fig. [4j shows that, for nrp=6, too many points have 
been cut. The respective profiles are not very different 
and a significant difference is only observed in the «e 
field, close to the impcllor. This small difference is thus 
able to change the onset about 15%. All these values are 
close to the LDV onset of R c m =42.5. 

In Fig. 10 we furthermore compare the isodensity sur- 



faces of the magnetic energy of two runs for LDV and 
SPIV data with identical parameters. SPIV calculations 
have been done on a much finer grid in axial direction. 
Both facings look very similar, which is confirmed by 
quantitative measurements. 



IV. KINEMATIC DYNAMO SIMULATIONS 



C. Test of the code and the preparation method 



A. Results 



Several kinds of tests have been performed in order 
to check the numerical code as well as the preparation 
methods of the SPIV data described in section IIII Al 

Ohmic diffusion. — The magnetohydrodynamic equa- 
tions ([2]) reduce for vanishing velocity field (u = 0) to 
the ohmic diffusion equation in cylindrical geometry 
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= Ab, 



(5) 



which has analytical eigenfunctions consisting of Bessel 
functions (T2]. We imposed these eigenfunctions in the 
code with a subsequent monitoring of the exponential 
decay of the amplitude. In all the investigated cases the 
theoretical eigenvalue has also been observed in the nu- 
merical run with high accuracy. 

Comparison with analytical flow results. — We used 
also the MND flow, an analytical velocity field studied 
in [T2] in order to ckeck our code. For the field given by: 



-r(l — r) 2 (l + 2r) cos(7rz), 



(6a) 



U0 



(l-r)(l + r-5r 2 )sin(7rz/2), (6b) 
4er(l — r) sin(7rz), (6c) 



with aspect ratio L = 2, for V = 0.776 and e = 0.743, the 
expected critical magnetic Reynolds number is 58, which 
we could confirm with our code. 

Comparison with LDV data. — In order to check our 
preparation procedure of the SPIV measurements we 



Onset of dynamo action. — We have investigated the 
configurations in which the two propellers are exactly 
counter-rotating, in (+) and in (-) direction in setups 
with and without annulus. The computations have been 
made with nrp=2 to 6, and the critical Reynolds number 
is found by computing the growth rate for different mag- 
netic reynolds numbers and the condition o-{R^ n ) = 0. 
The results are presented in Fig. |ll| and [T2] and summer- 
ized in Table H 

In the presence of an annulus in the shear layer, the re- 
sults show that the direction of rotation changes the crit- 
ical magnetic Reynolds number by a factor 3 : R^ ~ 36 
for the (+) case and R c m ~ 113 for the (-) case. There is 
an ever larger difference between the (+) and (-) rotation 
case when the annulus is absent: the dynamo threshold 
is R c m ~ 45 for the (+) case, and R c m > 500 in the (-) 
case. 

These differences can be traced back to a difference in 
the poloidal to toroidal ratio T, that is much lower in the 
(-) case (0.62 with annulus and 0.49 without annulus) 
than in the (+) case (0.8 with annulus or without annu- 
lus). This is in agreement with previous studies showing 
that the optimum for dynamo growth in von Karman 
type flows is around T sa 0.75 [T3]. Fig. 13 exhibits the 
influence of T on the magnetic growth rate a in the (-) 
case without annulus: the growth rate of the synthetic 
fields increases very sharply with T between T = 0.5 and 
T = 0.8. 

The insertion of the annulus exhibits two effects: (i) in 
the (+) case, it decreases slighly (15%) the critical mag- 
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netic Reynolds number while T remains unchanged; (ii) 
in the (-) case, it increases the poloidal to toroidal ratio 
r towards the range where dynamo action is available. 

In the (+) case, both mean poloidal and toroidal flow 
increase when the annulus is inserted, but their ratio T 
does not change within the error bars. However, with the 
annulus, the axial gradient of the toroidal flow is larger 
in the region of the shear layer (cf. Fig. [7|): the gradi- 
ent seems to be "less spread out" by the coherent struc- 
tures. In fact, the annulus has been intruduced in the 
VKS experiment to have an effect on the turbulent fluc- 
tuations, in particular on the large scale fluctuations, by 
stabilizing the shear layer vortices in the annulus plane. 
Indeed, the mean value of the normalized kinetic energy 
6(t) decreases from 5 = 2.02 to 1.48 and its variance 82 
decreases from 0.18 to 0.12 (cf. [H]). This means that 
in the configuration with the annulus, the instantaneous 
flow is much closer to the mean flow. This quantity has 
not been measured in the Riga and Karlsruhe experi- 
ments, but one can suspect that the instantaneous flow 
is very close to the mean flow, which explains the very 
good agreement between kinematic dynamo simulations 
and experiments for these configurations. This is not the 
case for the VKS experiment, in which turbulent fluctua- 
tions are present at all scales. The annulus thus appears 
as a mean to reduce this problem, in particular for the 
large scales. 

In the (-) case, T has a low value without annulus, be- 
cause both the toroidal part of the flow has increased and 
the poloidal part decreased. This low value is clearly out 
of the range for dynamo action in these flows. Moreover, 
the turbulent fluctuations are more intense in this sense 
of rotation and the corresponding values of 8 and 82 are 
larger than in the (+) case. The insertion of the annulus 
has two effects. First, it increases the poloidal to toroidal 
ratio, opening the possibility to get dynamo action, but 
still at large R c m . Second, it reduces 8 from 2.22 to 1.50 
and its variance 5% from 0.24 to 0.11 (cf. [H]), leading to 
values similar to the (+) case with annulus. 

Magnetic growing mode. — The fastest growing mode 
we observe is always am — 1 mode. Recall that since our 
configuration and flow are axisymmetric by construction 
and there are no external field source, the generation of 
a ttl — mode is forbidden by Cowling's theorem [31 j . 
Sections of the magnetic field obtained in the different 



configurations are presented in Fig. [T4jand Fig. 15 The 
configuration corresponding to (+) case without annu- 
lus is in agreement with the results of [T3]. The mag- 
netic field in the (+) case with annulus reveals a sharper 
field.... XXXXXXXX 

In the (-) case without annulus, the synthetic field 
with T = 0.8 looks very similar to the (+) case. On 
the contrary, the (-) case with annulus, that corresponds 
to a larger i?^ is very different from the 3 other cases.. 
XXXXXXXXXXXX 



B. Comparison with experiments 

Onset of dynamo action. — When comparing these 
kinematic simulations results with experimental one, one 
observes that if the exact values are note recovered, the 
order of magnitude of the dynamo threshold R c m as well 
as the general trend are correctly reproduced: with an- 
nulus, R c m is larger for the (-) case than for the (+) case, 
and there is no dynamo in the (-) case without annulus 
in the available experimental R m range. 

The fact that the thresholds obtained in simulations 
are different from the experimental ones is not surpris- 
ing: the simulations presented in this paper have been 
performed with real mean flows but in an over-simplified 
VKS configuration. First of all, the lateral boundary 
conditions (shell of sodium at rest) have been taken into 
account but not the end boundaries (sodium behind the 
disks). However, the effect is in general to increase i?^ 
f[20j). whereas the experimental values are smaller than 
the numerical ones. Then, the effect of the soft-iron im- 
pellers (magnetic permeability ~ 100) has also been dis- 
carded. Using the synthetic MND flow [12] to describe 
the von Karman mean axisymmetric flow, it has been 
shown that infinite permability boundary conditions are 
able to reduce the dynamo threshold by 5% to 10% |19j . 
Finally, and perhaps more importantly, turbulence is not 
adressed in the simulations, whereas the Reynolds num- 
ber of the VKS experiment is between 10 6 and 10 7 [5]. 
Only the mean axisymmetric flow is used in the simula- 
tions : the large scale vortices of the shear layer as well 
as the blades vortices are not considered, neither small 
scale turbulence. Non axisymmetric fluctuations, such as 
those generated by the flow ejected by the blades of the 
impellers, could e.g. generate an a effect and give rise to 
an alpha-omega type dynamo |15j . The addition in the 
induction equation of an empirical alpha term localized 
in the vicinity of the impellers displays an axial dipole 
but the value of a required to obtain the experimental 
threshold appears to be unrealistic [TO] . 

Very recent kinematic simulations with synthetic flows 
have been reported in two directions. First, by adding 
a modelization of the non axisymmetric vortices ejected 
by the blades of the propellers to the synthetic MND 
mean axisymmetric flow, |32j has evidenced dynamo ac- 
tion with an axial dipole or quadrupole for R c m ~ 34. 
Second, by adding a localized permeability distribution 
that mimics the shape of the impellers disk and blades, 
[33] has found an equatorial dipole for ~ 60 without 
a effect and an axial dipole for R c m ~ 32 when including 
an a effect with a realistic value. 

Magnetic growing mode. — In all the studied cases the 
structure of the kinematic neutral mode is a m = 1 equa- 
torial dipole mode. In the VKS experiment, the structure 
of the saturated has been reported in the (+) rotation 
sense, with annulus, and found to correspond mainly to 
a m = mode [6] [9]. It would be interesting to have 
detailed measurements of the saturated mode structure 
in the VKS experiment, for the different configurations 
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and in particular in the reverse rotation direction. One 
can e.g. wonder whether or not a m = 1 component is 
present besides the axial dipole. 

As reported above, recent simulations have exhibited 
an axial dipole (or quadrupole) with synthetic flows and 
particular boundary conditions, and it would be of great 
interest to study the effect of the annulus and the direc- 
tion of rotation on these results. 



V. CONCLUDING REMARKS 

Through kinematic simulations using SPIV velocity 
fields measured in a water model experiment, we have 
shown that the rotation direction of the non axisymmet- 
ric impellers in the VKS experiment has a great impact 
on the threshold of the dynamo generated by the mean 
flow only. By changing the rotation direction, one can 
increase the dynamo threshold by a factor 3 in the case 
with annulus (1.4 in the VKS experiment), and a fac- 
tor larger than 10 in the case without annulus (larger 
than 1.6 in the VKS experiment). Even if the exact ex- 
perimental results are not reproduced, the variation of 
the dynamo threshold with the direction sense is quali- 
tatively reproduced by the kinematic simulations. 

Since our kinematic simulations do not take into ac- 
count the magnetic permeability due to iron disks-an ef- 
fect that can decrease the threshold by a factor 2 (cf. 
[T^] ) or generate a m=0 mode with an a effect [33], we do 
not expect our results to be more than qualitative. Fur- 
thermore, our linear simulations do not take into account 



the back reaction of the Lorentz force, that can enforce 
emergence of the m=0 mode as a secondary bifurcation 
of the original m=l mode |17j . Finally, non axisymmet- 
ric turbulent fluctuations can play an important role in 
the generation of a magnetic dipole [H3 H2] . 

Nevertheless, we have shown that the increase (or the 
absence) of the dynamo threshold observed in given con- 
figurations of the VKS experiment can have two origins: 
(i) the poloidal to toroidal ratio T, a parameter previ- 
ously identified as crucial for dynamo action [14]. (ii) the 
normalized kinetic energy 8 which quantifies the distance 
of the instaneous flow to the mean flow. The parameter 
r is certainly important for the alpha-omega mechanism 
proposed in [15] . giving the relative importance of the a 
effect generated by the blades vortices (poloidal recircu- 
lation) vs. the omega effect due to differential rotation 
(toroidal flow). 

Detailed experimental study of the propellers vortices 
dynamics with the rotation direction would represent 
an important step to validate the numerical predictions. 
Work is in progress in this direction. 
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TABLE I: Comparison of experimental and numerical critical 
magnetic Reynolds numbers for the different VKS setups. An 
interval for R^sim is given when different nrp have been used 
(cf. Fig. |9]and 11 1, otherwise nrp=4. (*) denotes an artificial 
numerical field with a prepared F not available in experiments. 
The values of the mean normalized kinetic energy 8 and its 
variance 82 are from |llj 




FIG. 1: (a) Sketch of the simulated part of the VKS experi- 
ment (cross section) including both impellers (11,12) and an 
annulus (A) at mid height. The horizontal axis corresponds 
to the z axis. The container radius R is taken as length scale. 
The factor w denotes the thickness of the conducting layer 
(grey area) of sodium at rest in the VKS experiment. The 
layers of sodium between the impellers and the end walls of 
cylinder have not been taken into account in this work, (b) 
Top view of the TM 73 impellers used in the VKS experi- 
ment [14]. The (+) (resp. (-)) rotation sense corresponds to 
a counterclockwise (resp. clockwise) rotating impeller. 
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FIG. 2: Dimensionless velocity field measured by SPIV in 
VKE and symmetrized as well as spline extrapolated (nrp=4) 
for kinematic dynamo simulations (see text for details), (a) 
(+) rotation and no annulus (NA); (b) (+) rotation and pres- 
ence of an annulus (WA); (c) (-) rotation and no annulus 
(NA); (d) (-) rotation and presence of an annulus (WA); The 
coloring displays the velocity field perpendicular to the plane, 
whereas red (blue) denotes a positive (negative) sign in cylin- 
drical coordinates (r,Q,z). The white area for radii r > 1 
shows the conducting layer at rest (w — 0.6) with u = 0. 
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FIG. 3: Same as in Figure [2] (+) rotation has been used and 
there is no annulus. Part (a) indicates a section for fixed z = 
Z2 = —0.56, whereas in (b) an arbitrary cut of the rotational 
invariant field for fixed O is displayed. The arrows in (b) 
illustrate the section planes for the profiles shown in figures [4] 
and [5] 
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z z 

FIG. 4: Profiles of the velocity fields in cylindrical coordinates 
u = (u r , «e, u z) for fixed radii as explained in Fig. [3] A setup 
with (+) rotation has been used and there is no annulus. The 
left (right) column shows a section far away (nearby) the vor- 
tex center. Data obtained by processed LDV |Tj[] (circles) are 
compared with SPIV measurements (crosses). The lines show 
processed SPIV data used in the kinematic dynamo code for 
different spline extrapolation procedures, see text for details. 
All fields are in arbitrary units and scaled to the maximal am- 
plitude of the LDV data for comparison reasons. Note that in 
part (f ) the the scaling is expanded and the SPIV data values 
are very small. 
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FIG. 5: Profiles of the velocity fields in cylindrical coordinates 
u = (u r , U0,u z ) for fixed axial heights as explained in Fig. [3] 
and in caption of Fig. [4] (+) rotation has been used and there 
is no annulus. For these cuts (far away from the blades) there 
is no difference in the different spline extrapolations. Fields 
are shown in planes capturing their maximal value: (a) u r in 
the inflow area separating the vortex pair at z = Z\ = 0; (b) 
and (c) lie and u z in the vortex center at z — z% — —0.56. 
Data obtained by processed LDV [14] (circles) are compared 
with SPIV measurements (crosses) and processed data (line). 
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FIG. 6: Same as in Fig. [4] for a setup with (+) rotation and 
the presence of an annulus at mid height. 




FIG. 7: Comparison of the profiles of the velocity fields cor- 
responding to the 4 studied configurations for fixed radii as 
explained in Fig. [3] (+) rotation without annulus (contin- 
uous line); (+) rotation with annulus at mid height (dashed 
line); (-) rotation without annulus: artificial numerical field 
with F = 0.8 (dotted line); (-) rotation with annulus at mid 
height (dashed-dotted line) 




FIG. 8: Comparison of the profiles of the velocity fields corre- 
sponding to the 4 studied configurations for fixed axial height 
as explained in Fig. [3] (+) rotation without annulus (contin- 
uous line); (+) rotation with annulus at mid height (dashed 
line); (-) rotation without annulus:: artificial numerical field 
with r = 0.8 (dotted line); (-) rotation with annulus at mid 
height (dashed-dotted line) 




FIG. 9: Magnetic energy growth rate a vs. R m for a setup 
with a conducting layer (w = 0.4) and a flow obtained with 
(+) rotation and without annulus, F = 0.8. Full (open) sym- 
bols refer to SPIV (LDV) data. For SPIV the results for 
different spline fits are compared, cf. figures [4] and |5j respec- 
tively. 
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FIG. 10: Comparison of the isodensity surface of magnetic 
energy (50% of the maximum) for LDV (a) and SPIV (b). 
The figures show the neutral mode close to onset (R m = 44) 
with conducting layer (w = 0.6). A setup with (+) rotation 
and without annulus has been used. For SPIV the spline 
extrapolation with nrp=4 was used. The full line indicates 
the cylinder axis. Note, for x > (resp. x < 0) the angle 
is (resp. n). 




FIG. 11: Same as in Fig. [9]for a setup with (+) rotation and 
the presence of an annulus at mid height, V — 0.8. 
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FIG. 12: Magnetic energy growth rate a vs. R m for different 
setups with conducting layer (w = 0.4). WA (resp. NA) 
denotes configurations with (resp. without) annulus. For all 
cases nrp=4 has been used. Note that the case (-) rotation 
without annulus corresponds to an artificial numerical field 
with r = 0.8 not available in experiments. 
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FIG. 13: Influence of the poloidal to toroidal ratio Y on the 
magnetic growth rate a. T = 0.488 is the experimental value 
measured by SPIV in the experiment. The other cases cor- 
respond to artificial fields. A setup without annulus and (-) 
rotation has been used. For the spline extrapolation we fixed 
nrp=4. 
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FIG. 14: Sections of the magnetic field b for different exper- 
imental setups, here for fixed 6 = 0. All configurations are 
with conducting layer (w — 0.4), whereas WA (resp. NA) 
specifies setups with (resp. without) annulus. Arrows cor- 
respond to components lying in the cut plane and the color 
code to the component transverse to the cut plane, where red 
(resp. blue) denotes a positive (resp. negative) sign. Param- 
eters are (a) R m = 44, Y = 0.8, (b) R m = 40, F = 0.8, (c) 
R m = 60, T = 0.8 (artificial field), (d) R m = 115, F = 0.62. 
For all cases nrp=4 has been used. 




FIG. 15: Same as in Fig. [TT] for 9 = tt/2 



